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A laser, be it an optical laser or an atom laser, is an open quantum system that produces a 
coherent beam of bosons (photons or atoms respectively). Far above threshold, the stationary state 
pan of the laser mode is a mixture of coherent field states with random phase, or, equivalently, a 
Poissonian mixture of number states. This paper answers the question: can descriptions such as 
these, of p^s as a stationary ensemble of pure states, be physically realized? Here physical realization 
is as defined previously by us [H.M. Wiseman and J. A. Vaccaro, Phys. Lett. A 250, 241 (1998)]: 
an ensemble of pure states for a particular system can be physically realized if, without changing 
the dynamics of the system, an experimenter can (in principle) know at any time that the system is 
in one of the pure-state members of the ensemble. Such knowledge can be obtained by monitoring 
the baths to which the system is coupled, provided that coupling is describable by a Markovian 
master equation. Using a family of master equations for the (atom) laser, we solve for the physically 
realizable (PR) ensembles. We find that for any finite self-energy x of the bosons in the laser 
mode, the coherent state ensemble is not PR; the closest one can come to it is an ensemble of 
squeezed states. This is particularly relevant for atom lasers, where the self-energy arising from 
elastic collisions is expected to be large. By contrast, the number state ensemble is always PR. As 
the self-energy x increases, the states in the PR ensemble closest to the coherent state ensemble 
become increasingly squeezed. Nevertheless, there are values of x for which states with well-defined 
coherent amplitudes are PR, even though the atom laser is not coherent (in the sense of having a 
Bose-degenerate output). We discuss the physical significance of this anomaly in terms of conditional 
coherence (and hence conditional Bose degeneracy). 
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I. INTRODUCTION 

In elementary presentations of quantum optics it is 
more or less an axiom that a laser field is represented by a 
coherent state \a). Recently, it has been argued that this 
representation is a fiction, albeit a convenient one [y. 
The essential argument is that no commonly employed 
process at optical frequencies produces an electric field 
having a non-zero average amplitude. While this point 
of view is certainly defensible [g[ , it perhaps obscures the 
fact that there is something special about laser light. 

In Ref . 0] , one of us argued that what is special about 
laser light is that it is well approximated by a noiseless 
classical electromagnetic wave. Four quantitative crite- 
ria were given, none of which require a mean field, so 
there is no dispute with Ref. [Q. The least familiar, and 
so most important, of these criteria is that the output 
flux of the laser (bosons per unit time) must be much 
greater than its spectral linewidth. Put another way, the 



coherence time of a true laser must be much greater than 
the mean temporal separation of photons in the output 
beam. This is typically satisfied by many orders of mag- 
nitude in optical lasers, but is not satisfied by ordinary 
thermal sources. 

This concept of quantum coherence is quite distinct 
from the elementary idea that a laser is in a coherent 
state. Indeed, it is compatible with theoretical models for 
typical laser processes Mpl , which imply that the state of 
the cavity mode for a laser far above threshold is a mix- 
ture of coherent states of all phases. That is to say, the 
stationary state matrix of the laser mode can be written 
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\aW^){\aW^\. 



(1.1) 



where \a\^ = ^ is the mean number of photons in the 
laser. 



It would be tempting to interpret Eq. (1.1) to mean 
that the laser really is in a coherent state I |a|e*'^) of def- 
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inite phase (/), but we don't know what that phase is. 
However, this temptation must be resisted because the 
stationary state matrix can also be written 



^e-'^^|n)(n| 

n=0 



(1.2) 



which would seem to imply that the laser really is in a 
number state \n), but we don't know which number it is. 

The "unknown coherent state" description and the 
"unknown number state" description are mathematically 
equivalent representations of the stationary state matrix 
Pss- However, in the physical context that pss is the sta- 
tionary state of an open quantum system in dynamical 
equilibrium, the two representations are not physically 
equivalent. This idea is at the heart of this paper and the 
following paper |g] . In this paper we investigate whether 
these, and other pure state ensembles are physically real- 
izable. We will show that under some circumstances, the 
"unknown coherent state" description is not physically 
realizable, in contrast to the "unknown number state" 
description, which is. In the following paper we look at 
the question of the how robust the ensembles are. We find 
that even among physically realizable ensembles, a phys- 
ical distinction may be drawn based upon the survival 
time, the average time that a member of the ensemble 
remains close to its original state when left to evolve un- 
der the system dynamics. Both of these concepts, the 
physical realizability of pure-state ensembles, and the ro- 
bustness of such ensembles, were introduced in an earlier 
work by us |Q . 

Before proceeding further, it is necessary to clarify 
what we mean by "physically realizable" (PR). A sta- 
tionary pure state ensemble of a given system is PR if 
it is possible, without altering the dynamics of the sys- 
tem, to know that its state at equilibrium is definitely 
one of the pure states in the ensemble. Of course which 
pure state cannot be predicted beforehand. It may seem 
contradictory to say that the system at equilibrium is 
mixed, but that nevertheless we can know it to be in a 
pure state. The resolution is that, by monitoring the sys- 
tem's environment, the system state can, under suitable 
circumstances, be collapsed over time into a pure state. 
Being simply an example of a quantum measurement, 
this process, called an unraveling [p|, will be stochas- 
tic. On average, the system evolution is not changed and 
the ensemble of pure states produced by the unraveling 
is guaranteed to be equivalent to the equilibrium mixed 
state. 

From this description it should be apparent that the 
question of whether an ensemble is PR or not can- 
not be determined from the stationary mixed state pss- 
Rather, it depends upon the dynamics (reversible and ir- 
reversible) that produced the stationary state. Indeed, 
the unraveling to a pure state is realized by monitoring 
the environment of the system, the same environment 
that produces the irreversible dynamics of the system. It 
would not be justifiable to introduce some new reservoir 



to allow a new measurement to be made. Even if that 
did not change the stationary state of the system (such 
as would be the case for adding a QND-boson number 
measurement to a laser), it would change the dynamics 
of the system, and hence one would be investigating a 
different system. 

The fact that different dynamics can lead to the same 
stationary mixed state is easy to see for the case of a 
laser. Any process that commutes with boson number 
will not alter the stationary laser state pss, since its eigen- 
states are the number states, as shown by Eq. (1.2). An 
example of an irreversible process that commutes with 
boson number is phase diffusion. This is relevant to all 
current lasers, which have some phase diffusion in excess 
of the standard limit (although see Ref. ^J for theoret- 
ical proposals for lasers that have phase diffusion below 
the standard limit). There are also reversible processes 
that commute with boson number, such as degenerate 
four-wave mixing. While this dynamics is unimportant 
in most optical lasers, it is expected to be very significant 
in atom lasers. 

An atom laser is a device that produces an output 
beam of bosonic atoms analogous to an optical laser's 
beam of photons [B. The idea for an atom laser was 
published independently by a number of authors [p^-p^ , 
shortly after the first achievement of Bose-Einstein con- 
densation (BEG) of a dilute atomic gas p4|-p6[. There 
have since been some important experimental advances 
in the coherent release of pulses |l^,^ql and beams [19 20 
of atoms from a condensate. Because the condensate is 
not replenished in these experiments, the output coupling 
cannot continue indefinitely, so these devices are only the 
first steps towards achieving a CW atom laser. 

Even though the atoms in the current BEG experi- 
ments are weakly interacting in the sense of forming a gas 
rather than a liquid, elastic collisions may dominate the 
dynamics of the condensate. This self-interaction does 
not directly alter the number of atoms in the conden- 
sate, and is analogous to four-wave mixing, (that is, a 
X'^' nonlinearity) , in optics. In this paper we show that 
the presence of this nonlinearity has an enormous influ- 
ence on what ensembles of pure states are physically re- 
alizable. It also determines the laser linewidth, and in 
this paper, we explore the connection, between the pres- 
ence of a PR coherent amplitude, and the coherence of 
the laser output. 

This paper is organized as follows. In Sec. II we ex- 
plain in detail our concept of physically realizable pure- 
state ensembles for open quantum systems. In Sec. Ill we 
present our atom laser model, including self-interactions 
and phase diffusion. In Sec. IV we apply the formalism 
of Sec. II to the atom laser model and set up the frame- 
work for calculating the PR ensembles. We calculate the 
PR ensembles in Sec. V and derive various scaling laws 
for the ensembles as a function of the self interaction and 
phase diffusion. We conclude in Sec. VI with a summary 
and a discussion of the interpretation and implications of 
our work. 



II. PHYSICALLY REALIZABLE ENSEMBLES 
A. The Master Equation 

Open quantum systems generally become entangled 
with their environment, and this causes their state to 
become mixed. In many cases, the system will reach an 
equilibrium mixed state in the long time limit. A CW 
laser or atom laser is a system of this sort, and we will 
restrict our consideration to such systems. It is common 
to refer to the environment of these systems as a reser- 
voir and, accordingly, we use both terms (environment 
and reservoir) interchangeably here. 

If the system is weakly coupled to the environmental 
reservoir, and many modes of the reservoir are roughly 
equally affected by the system, then one can make the 
Born and Markov approximations in describing the ef- 
fect of the environment on the system pi| . Tracing over 
(that is, ignoring) the state of the environment leads to 
a Markovian evolution equation for the state matrix p of 
the system, known as a quantum master equation. The 
most general form of the quantum master equation that 
is logically valid is the Lindblad form p2[ 



K 



P 



-i[H,p] + Y,'D[ck]p = Cp, 



fe=i 



where for arbitrary operators A and B, 
V[A]B = ABAl - {A^A, B}/2. 



(2.1) 



(2.2) 



If the master equation has a unique stationary state 
(as we will assume it does), then that is defined by 



Cps: 



0. 



(2.3) 



This assumption requires that C be time- independent. In 
many quantum optical situations, one is only interested 
in the dynamics in the interaction picture, in which the 
free evolution at optical frequencies is removed from the 
state matrix. Indeed, for quantum systems driven by a 
classical field, it may be necessary to move into such an 
interaction picture in order to obtain a time-independent 
Liouvillian superoperator C 

The stationary state matrix pss can be expressed as an 
ensemble of pure states as follows: 



Ps: 



E 



^71-^71-) 



(2.4) 



where the P„ are rank-one projection operators 

P„-|^n)(V'n|, (2.5) 

and the p„ are positive weights summing to unity. The 
(possibly infinite) set of ordered pairs. 



S = {(P„,p„):n=l,2,...}, 



(2.6) 



we will call an ensemble E of pure states. Note that 
there is no restriction that the projectors P„ be mutu- 
ally orthogonal. This means that there are continuously 
infinitely many ensembles E that represent pss- As noted 
in the introduction, only some of these are physically re- 
alizable. 



B. Unravelings 

In the situation where a Markovian master equation 
can be derived, it is possible (in principle) to continually 
measure the state of the environment on a time scale large 
compared to the reservoir correlation time but small com- 
pared to the response time of the system. This effectively 
continuous measurement is what we will call "monitor- 
ing" . In such systems, monitoring the environment does 
not disrupt the system-reservoir coupling and the system 
will continue to evolve according to the master equation 
if one ignores the results of the monitoring. 

By contrast, if one does take note of the results of mon- 
itoring the environment, then the system will no longer 
obey the master equation (except on average). Because 
the system-reservoir coupling causes the reservoir to be- 
come entangled with the system, measuring the former's 
state yields information about the latter's state. This 
will tend to undo the increase in the mixedncss of the 
system's state caused by the coupling. 

If one is able to make perfect rank-one projective (i.e. 
von Neumann p3|) measurements of the reservoir state, 
with negligible time delay from when it interacted with 
the system, then the system state will usually be col- 
lapsed towards a pure state. However this is not a process 
that itself can be described by projective measurements 
on the system, because the system is not being directly 
measured. Rather, the monitoring of the environment 
leads to a gradual (on average) decrease in the system's 
entropy. 

If the system is initially in a pure state then, under 
perfect monitoring of its environment, it will remain in 
a pure state. Then the effect of the monitoring is to 
cause the system to change its pure state in a stochastic 
and (in general) nonlinear way. Such evolution has been 
called a quantum trajectory pi, and can be described 
by a nonlinear stochastic Schrodinger equation [P4|-p6[ . 
The nonlinearity and stochasticity are present because 
they are a fundamental part of measurement in quantum 
mechanics. 

Although a stochastic Schrodinger equation is concep- 
tually the simplest way to define a quantum trajectory, 
in this work we will instead use the stochastic master 
equation (SME) ||-||. 

This has four general advantages. First, it can de- 
scribe the purification of an initially mixed state. Sec- 
ond, it can easily be generalized to describe the situa- 
tion where not all baths are monitored perfectly, and the 
conditioned state never becomes pure (as we will con- 



sider in Sec. VI). Third, it is easier to see the relation 
between the quantum trajectories and the master equa- 
tion that the system still obeys on average. Fourth, the 
form of the SME is invariant under stochastic U{1) trans- 
formations of the state vector, which can radically alter 
the appearance (but not the substance) of the stochastic 
Schrodinger equation Q. 

Assuming that the initial state of the system is pure, 
the quantum trajectory for its projector P will be de- 
scribed by the SME 



dP = dt[C + Uit)] P. 



(2.7) 



Here C is the Liouvillian superoperator from the master 
equation, and U is a stochastic superoperator which is, 
in general, nonlinear in its operation on P. It also de- 
pends on the operators Ck as defined in Eq. (2.1), and is 



constrained by the following two equations, which must 
hold for arbitrary rank-one projectors P 



{P,{C+U)P} 



dt[UP][UP] = (C- 
E[UP] = 0. 



-U)P, 



(2.8) 
(2.9) 



The first of these properties ensures that P + dP is 
also a rank-one projector; that is, that the state remains 
pure. The second ensures that 



dE[P] = CE[P]dt, 



(2.10) 



where E denotes the ensemble-average with respect to the 
stochasticity oiU. This stochasticity is evidenced by th e 
necessity of retaining the term dt[UP][UP] in Eq. (2 



Because the ensemble average of the system still obeys 
the master equation, the stochastic master equation (or 
equivalently the stochastic Schrodinger equation) is said 
to unravel the master equation |q] . It is now well-known 
P4| that there are many (in fact continuously many) 
different unravelings for a given master equation, cor- 
responding to different ways of monitoring the environ- 
ment. 

For simplicity we will call U an unraveling. Each un- 
raveling gives rise to an ensemble of pure states 



E' 



= {(P,'/,p^):n = l,2, 



C. Continuous Markovian Unravelings 

To determine whether an ensemble £^ is a PR ensemble 
E'^ requires a search through the set, call it J, of all possi- 
ble unravelings U. This set is extremely large. Although 
the stochasticity in the superoperators U can always be 
written in terms of quantum jumps, these jumps range in 
size from being infinitesimal, to being so large that the 
system state after the jump is always orthogonal to that 
before the jump [p5| . 

Another complication is that the unraveling need not 
be Markovian, even though the master equation is. It 
might be thought that the measurement must be Marko- 
vian since it must obtain full information from the en- 
vironment immediately after it has interacted with the 
system in order that the conditioned system state re- 
main pure. This rules out spectral detection, for exam- 
ple, where the conditioned system state is not pure be- 
cause it is entangled with the state of the spectral filters 
|36| . However, the way in which the measurement obtains 
information from the environment may not be indepen- 
dent of the past history of the system. For example, 
the parameters defining the measurement may depend 
on previous measurement results, resulting in an adap- 
tive measurement, as discussed in Ref. p7[|. 

From these considerations we see that a search over all 
possible unravelings would not be practical. Thus it is 
useful, to consider a smaller (but still continuously infi- 
nite) set D containing only unravelings that are continu- 
ous and Markovian (CM). A continuous (but not diffcr- 
entiable) time evolution arises from infinitely small (and 
infinitely frequent) jumps |^,0. In this case the proba- 
bility distribution for the pure states obeying the SME 
satisfies a Fokker-Planck equation. On this basis it has 
been argued that these unravelings are the natural ones 
to consider for quantum systems expected to show quasi- 
classical behavior |33] . The measurement will be Marko- 
vian provided the measurement parameters Ujk (defined 
below) are constants. 



For the general master equation (2.1) the elements U 
of D can be written as RBSl 



K 



•}, 



(2.11) 



u{t)dt^j2'^\-(^^k(t) 



Ck 



(2.13) 



fc=i 



where P^ are the possible pure states of the system at 
steady state, and p^f are their weights. For master equa- 
tions with a unique stationary state pss, the SME (2.7) 



is ergodic over E^ and pjf is equal to the proportion of 
time the system spends in state P^ . The ensemble E^ 
represents pss in that 



Vp^F^=p,„ 



(2.12) 



as guaranteed by Eq. ( 2.10| ) 



Here Ti-IA] is a nonlinear superoperator defined, for arbi- 
trary operators A and i?, by 



n[A]B = AB + BA^ - Tr[AB + BA'^]B, 



(2.14) 



and the dWk{t) are the infinitesimal increments of a com- 
plex multi-dimensional Wiener process [ p8[ satisfying 



E[dH/fe] = 
dWj{t)dW^{t) =dt6jk, 
dWj{t)dWkit) =dtujk. 



(2.15) 
(2.16) 
(2.17) 



The only condition on the complex numbers Ujk — Ukj 
is that the corresponding complex symmetric matrix u 
must satisfy [BSl 



lull < 1. 



(2.18) 



This conies from the requirement that the following ex- 
pression must be non-negative 



E 



^ ZkdWk + c.c I 
^ k ) 



(2.19) 



for an arbitrary iiT- vector of complex numbers z. 

Some insight into the measurement parameters Uj^ 
may be found by considering the simple case with one 
irreversible term; that is, K ~ \ so that there is just one 
complex number u in Eq. (2.17). For specificity, say the 
system is an optical cavity with annihilation operator a, 
damped through one end mirror with decay rate k. Then 
the continuous Markovian unravelings correspond to two 
independent homodyne detection apparatuses H , each of 
efficiency 1/2. If the local oscillator phases are Q\ and Qi 
then u = (e^'^i -I- e'^'-'^'^)/2. The two photocurrents hit) 
and I2 (t) , normalized to have unit shot noise, are given 

by [III 

Ip{t)dt^ ^/^{e-''^''a + e'^''a'')dt + dWp{t), (2.20) 

where dWi and dW2 are independent Wiener increments. 
We can combine the photocurrents to make a complex 
signal 

J{t)dt = [e*^i/i(t)di -I- e'^''l2{t)dt]/V2 (2.21) 

^ y^{a + ua^)dt + dW{t), (2.22) 

where dW{t) = [e''^^dWi{t) + e''^^dW2it)]/V2 is a com- 
plex Wiener increment satisfying 

dW*(t)dW{t)^dt, dW{t)dW{t)=udt. (2.23) 

That is, it has the same correlations as the dW{t) occur- 
ring in the stochastic master equation, and is in fact the 
same noise process. 

If the two local oscillator phases are chosen to be iden- 
tical then |u| = 1 and both apparatuses measure the same 
quadrature of the cavity mode. If they are chosen to be 
in quadrature, with Oi — 62 = 7r/2, then u = Q and two 
orthogonal quadratures are measured each with efficiency 
1/2. In general < |w| < 1, and for any u 7^ 0, differ- 
ent amounts of information are obtained about the two 
cavity-field quadratures. The information gained tends 
to reduce the cavity field to a state with correspondingly 
different quadrature uncertainties. This gives an idea as 
to how different unravelings can give rise to different en- 
sembles. 

For a master equation with K Lindblad terms the 
problem of finding the ensembles that are physically real- 
izable by some continuous Markovian unraveling (CMU) 



reduces to determining the boundary {ujk : |lu|l = 0} of 
a region in K{K -\- l)-dimensional Euclidean space. Even 
for a moderately sized K (for example K — i \s needed 
for the atom laser problem), this is a surprisingly large 
space, which is difficult to search efficiently. For that 
reason we adopt in this paper a different search strategy, 
which will be explained in Sec. IV C . 



D. Quantum State Diffusion 

There is an interesting continuous Markovian unravel- 
ing, which has some special properties, for the case where 
Uij = |^,Q. In this case each complex Wiener pro- 
cess dW can be decomposed into real Wiener processes 
dW.dW^ as 



dWk = {dW^ -f idWl)/V2 



(2.24) 



such that dW^dW°^ 



dW^dW^ 



dWldW}> 



Sjkdt, and 
0. This unraveling is invariant under the 
complete set of linear transformations of the Lindblad 
operators, 



f/u 



(2.25) 



that leaves the master equation invariant. Here [/^^ is 
an arbitrary unitary matrix. 

This unraveling was introduced by Gisin and Percival 
| pO| , under the name of quantum state diffusion (QSD), 
as a microscopic model for decoherence. In the opti- 
cal context, it has been interpreted as the unraveling 
resulting from heterodyne detection |g^ or from equal- 
efficiency homodyne detection of orthogonal quadratures 
(as discussed above), although it can also arise in atomic 
detection schemes as well [^. It has been suggested by 
Rigo and Gisin |35| that the QSD unraveling is a natural 
way to discover the classical limit for aquantum system. 
Along similar lines, Diosi and Kiefer ETl have suggested 
that the QSD unraveling is the most robust unraveling, 
or close to the most robust unraveling (see the following 
paper Q for a detailed discussion of this concept). Thus, 
as well as considering the set of all ensembles physically 
realizable from CMUs, we will also pay particular atten- 
tion to the ensemble arising from the special instance of 
QSD. 



E. Discontinuous Unravelings 

Although most of our calculations are restricted to 
CMUs, there will be one occasion where we need to con- 
sider the following discontinuous unravelings of the mas- 
ter equation (2.1): 



U = J2^k, 



(2.26) 



where 



Uk{t)pdt ^ [dNk(t) ~ dtXkip)] 






Here the dNk{t) are pomt processes defined by 

dNj{t)dNk{t) ^ SjkdNkit) 

E[dNk] = Xk{p)dt = TT[pclck]dt 



(2.27) 



(2.28) 
(2.29) 



It is easy to verify that t his u nrave ling satisfies the nec- 



essary conditions of Eqs. (2.8) and (2.9). This unravehng 



(quantum jumps) is the most commonly used for numer- 
ical simulation of master equations |42] . 



III. THE (ATOM) LASER 

The system we wish to consider in this paper is the 
(atom) laser. As noted in the introduction, we take a 
laser to be a device that produces a coherent output, in 
the sense explained in Ref. [pj. An atom laser is thus a 
device that produces a coherent beam of bosonic atoms, 
analogous to the coherent beam of photons from an op- 
tical laser. 



A. The Master Equation 

A generic model for a laser was derived in Ref. [0. 
It describes a single-mode field having annihilation op- 
erator a, evolving under linear damping and nonlinear 
amplification. The nonlinearity in the amplification is 
due to depletion of the source (the gain medium in opti- 
cal lasers) and is essential for a coherent output to form. 
In the interaction picture, and measuring time in units 
of the decay rate, the master equation is 



P 



pV[a1] {A[< 



p + V[a]p. 



(3.1) 



The two terms on the right describe saturated gain and 
the decay due to the coupling of the laser mode to the 
output beam, respectively. Here Ug is the saturation bo- 
son number, /i is a (typically) large parameter, 2? is as 
defined in Eq. (|2.2|) and for arbitrary operators A and B, 



(3.2) 



A[A]B = [A^AB + BA^A]/2. 



For simplicity we take the limit where Us can be ig- 
nored compared to aa^ . Strictly this requires the limit 
Us ^ I, because the smallest eigenvalue of aa^ is 1. How- 
ever, for a laser at steady state the mean boson number 
is typically much greater than 1, and only boson num- 
bers close to the mean are occupied with any significant 
probability. In the above model the mean number is ap- 
proximately p — Hg in the limit of large u. H ence in the 
limit /i ^ rig, 1 we can ignore Us in Eq. (3.1). The resul- 
tant far-above-threshold laser master equation was first 
exphcated in Ref. [H. 



Having made this simplification we now introduce 
more terms into Eq. (3.1) in order to create a more real- 
istic model. First, we introduce a term describing phase 
diffusion. This will be present in optical lasers for all 
sorts of technical reasons such as thermal motion of the 
cavity mirrors. In an atom laser it may also be present 
for more fundamental reasons, such as collisions between 
uncondensed atoms (in the source modes) and atoms in 
the laser mode condensate. Treating this phase diffusion 
as a Markovian process, it is described by a Lindblad su- 
peroperator of the form NVla^a], where N is the phase 
diffusion rate in units of the decay rate. 

The second new term we introduce is peculiar to atom 
lasers: the self-energy of atoms in the condensate. This 
is described by a Hamiltonian equal to hG{a^af' , with 



C 



2'Kha^ 



d3r|V'(r)|^ 



(3.3) 



where ^(r) is the wavefunction for the condensate mode, 
tts is the s-wave scattering length, and k is the unit- 
valued decay rate of the condensate. Like the extra phase 
diffusion term, this term has no effect on boson number; 
it only affects the phase of the field. However it is strictly 
not a phase diffusion term, but rather a dispersive term. 
It would arise in an optical laser in a medium with a 
nonlinear refractive index. 

Putting the four terms (gain, loss, phase diffusion and 
self-energy) together, the total master equation is 



P 



= {pV[a)]A[c 



ti 



V[a\ + NV[a)a]) p 



iC[{a^a)\p]. 



(3.4) 



That this is of the Lindblad form follows from the iden- 
tity 



tl ii[„ti-i 



V[a?]A[a 



dqVla'e 



t -(jaat/Sl 



(3.5) 



The stationary solution is a Poissonian mixture of num- 
ber states with mean /i, just as expressed in Eqs. (1.1) 
and (|l|): 



Pss = / ^ IVpe^'^) (Vpe^^l = £ e-^^ \n) {n\ . (3.6) 



n=0 






B. The Linearized Master Equation 



The master equation ( |3.4[ ) is rather difficult to deal 
with because of the nonlinearities in both the gain term 
and the self-energy term. To make it more tractable we 
linearize this equation for a state localized about a mean 
field (a) = y//I. We make the replacement 



a = V7^+ {x + iy)/2 



(3.7) 



and get, to second order in x and y, 

p = (1/4) {D[x + ly] + (1 + v)V[x\ + V[y] 

+ n[i{xy + yx)l2^ixx']]p, (3.8) 



where 



u = ANfi > 0, X = 4/iC 



(3.9) 



and Ti. is the superoperator defined in Eq. (2.14), which 
here is serving as a convenient way to describe the Haniil- 
tonian evolution. We have ignored a contribution to the 
hnearized Haniiltonian that is proportional to a'^a as this 
simply indicates a frequency shift that can be removed 
in the interaction picture. 

To solve this master equation, we use the Wigner rep- 
resentation W{x,y) p^l. We make a Gaussian ansatz 



W{x,y) = exp 



M20M02 



M20M02 — Ml 



JX- fiiof 
2^20 



+ 



Pii{x - nio){y - uni) (y-Moi)' 



M20M02 



2^0 



27rW/i2oMo2 - Mil 



Substituting this into Eq. (3^) yields the following ODEs 
for the moments 



MlO — ^MlO: 

Moi = -XMio, 
M20 = -2^20 + 2, 
Mil = -Mil -XM20, 
Mb2 = -2xMii + 2 + iy. 

The solution is easy to find 

Mio(i) = Mio(0)w, 

Moi(i) = Moi(O) - XMio(0)(l - w), 

M2o(t) = M2o(0)w^ + 1 - w^, 

Mii(i) = Mii(0)w - x{l + u'[M2o(0) - 2] 

+ w''[l"M2o(0)]}, 
M02(t) = ^02(0) + (2 + i/)t - 2xMii(0)(l - w) 
+ 2x'{t+[M2o(0)-2](l-ii;) 

Here we are using the abbreviation w = e~* . 



(3.11) 


(3.12) 


(3.13) 


(3.14) 


(3.15) 


(3.16) 


(3.17) 


(3.18) 



(3.19) 



(3.20) 



C. Coherence 

Having solved for the dynamics of our (atom) laser 
model, we can now answer the question, is it a true laser? 
That is, does it satisfy the criteria for a coherent output 
as detailed in Ref. 0] . The first two criteria will be satis- 
fied provided the output coupling is realized in a suitable 



way. The next two relate to the quantum noise of the 
state, and depend upon the dynamics. 

First, the laser intensity should be well-defined. Al- 
though this criterion is strictly defined in terms of the 
output of the laser, it will be satisfied if the boson num- 
ber of the laser mode itself is well-defined. In the present 
case this is clearly so provided the mean number satisfies 



M>1, 



(3.21) 



as the ratio of the standard deviation to the mean is equal 
to 1/^. 

Second, the laser phase should be well-defined in the 
sense that the phase should stay approximately constant 
over the time between the emission of one boson and the 
next. With a unit damping rate, this time is equal to 
M~^. Rigorously, we require that the magnitude of the 
first order coherence function 



g(i)(t) = (at(i)a(0))/(ata) 



(3.22) 



remain close to unity for i = /i ^ . For the current system 
we can rewrite this expression as 



(3-10) g«(i)=M-^Tr[ate/*(ap3s)] 



(3.23) 
= p-^ f d^Ti [a^e^\a \y^e"^) (y^e'^^l)] (3.24) 

Now because £ is a phase-independent superoperator, 
the trace here is independent of (p. Thus the integral can 
be dropped and we can rewrite this as 



gW(i) = (l/a*)Tr[a 



te^* la) 



H] 



(3.25) 



= fi. Thus, the requirement that g^^'{t) ~ 1 
^ is exactly equivalent to requiring that the 



where \a\' 

for t — fi 

system, initially in a coherent state of mean number /z, 

still have a phase variance much less than unity after a 

time t — M^^- 

Without loss of generality we can take the initial coher- 
ent state to be jy^). Then Mio(O) = Moi(O) = Mii(O) = 
0, M20 = M02 — Ij and y is the phase quadrature. Assum- 
ing that the phase uncertainty remains relatively small, 
we can make the approximation 



y 
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From Eq. (3.17), the mean phase remains zero 



while the phase variance increases as 



(3.26) 



(3.27) 



4m 
Substituting t = /i^^ < 1 into Eq. ( ^.2C| ) yields 



(3.28) 



\t)) = 



l + (2 + t/)/i-i+xV'^ 



(3.29) 



For the phase to remain well-defined we require this to 
be much less than unity. Since we already require fi^ I, 
this gives the extra conditions 



3/2 






(3.30) 
(3.31) 



In a typical optical laser (and certainly in some models 
of atom lasers |1^] ) , i' ^ 1 . This means that excess phase 
diffusion dominates the intrinsic phase diffusion (which 
gives the 2 in the 2 + v term). In a typical atom laser, 
it is also likely that excess phase diffusion will dominate. 
However, as long as i^ <C /i^ the laser will remain coher- 
ent. Since v — 'iNfi, this is equivalent to the condition 



A^<M- 



(3.32) 



This expression places an upper bound on the phase dif- 
fusion rate N for the device to be considered a laser. 

For an optical laser any nonlinear refractive index is 
usually small and x ^ 1- For an atom laser x is lik ely to 
be much greater than one, as we will discuss in Sec. VI D. 
To be a true atom laser it is necessary for it to remain 
much less than /i^'^. Since x = ^^-C the phase coher- 
ence condition pla ces a n upper bound on the condensate 
self-energy in Eq. (3.3) of 



C< /x 



1/2 



(3.33) 



IV. UNRAVELING THE (ATOM) LASER 

We now wish to consider monitoring the environment 
of the laser in order to realize physically an ensemble of 
pure states. This would be very to difficult to do exper- 
imentally, as it would require monitoring all reservoirs 
for the device, including the source of bosons (the gain 
medium) and the sources of phase diffusion as well as the 
laser output. However in principal these things can be 
done providing the laser evolution is well-approximated 
by a Markovian master equation. 



A. Realizing the Number State Ensemble 

Before turning to continuous Markovian unravelings, 
we consider a discontinuous unraveling to show how the 



ensemble consisting of number states can always be re- 
alized. Using the atom laser master equation (3^) in 
the Lind blad form ( |3.5| ), we can apply the unraveling of 
Sec [IE, where the Lindblad operators are 



Co = a 

Cm — y Na*a, 



(4.1) 
(4.2) 



plus a continuum of Lindblad operators 

Cq = T^a^e-'?""*/^, for q G [0,oo). (4.3) 

Each of these operators either leaves the a number state 
unchanged, or turns it into another number state. Since 
the Hamiltonian Ca^ a^ aa also leaves a number state un- 
changed, it follows that if the system is initially in a 
number state, it will simply jump between number states 
under this unraveling. Moreover, it can be shown that 
an arbitrary initial state will tend towards some number 
state under this unraveling. In this way it is clear that 
the number state ensemble (1.2) can always be physically 
realized. 

B. The Continuous Markovian Unravelings 



As mentioned in Sec. |IQ, we are principally concerned 
with continuous Markovian unravelings. In this case, 
from the master equation (3.4), the SME is 



+ VJI dqn[dW*{t)a^ 
Jo 



,tg-9QQV2lp 



+ H[dWo{t)a]P + VNn[dWM{t)a^a]P 
-idt[C{a^af,P]. (4.4) 

Here dWo is a zero-mean white noise term. If we define 
Co(i) = dWoit)/dt we have 



nCo{t)(:o{t')]^S{t-t'). 



(4.5) 



and likewise for Cat and (q for each q. We say that these 
white noise terms are distinct because the cross terms are 
zero, for example 

nCoit)CN{t')]=0. (4.6) 

Now we wish to linearize. First note that 



Vm / dqCit)a^ exp(-gaat/2) ~ ^i / dqC{t)e-f^^^^^[l + {x + ly - iixq)/2^ 



c-numbcr+ I / dqiQ{t)c-'^''/^ + | / dqQ{t)e-^^'^'\l - q) 



= c-numbcr + [yQ{t) + xQ{t)] /2, 



(4.7) 

(4.8) 
(4.9) 



where C2(i) and CaC^) are disti nct c omplex normalized white noise terms as usual. 
Using this, we linearize Eq. ( [4.4|) as 



dP = (l/4)di {V[x + iy]p + (1 + '^)V[x]p 
+ (1/2) {n[dW^{t){x + ly)] + VT+ 



where we have defined 



a new white noise source 
\/iydWN{t). We could have 



^/TTl^dWiit) = dWait) 
obtained this result direc tly from the linearized form of 
the master equation (|3.8| ), but this derivation makes the 
physical origin of the noise terms apparent. 

The three complex white noise sources dWj = Q dt are 
distinct in the above sense that 

nQm{t')]^s,At-t')- (4-11) 

However they can still be correlated in the sense that 

EUt)Q{t')]=^vSit~t'), (4.12) 



f V[y]p + n[i{xy + yx)l2] + H[-ixx'']) P 

'^n[dw*{t)x] + n[dw;{t)y\] P. 



(4.10) 



where the Uij are constrained only by Eq. ( |2.1^ ). The 
(5-function in time in Eq. (4.12) is not required to re- 
produce the master equation. It is a consequence of our 
restriction to Markovian unravelings. 

Now it is a remarkable fact about the stochastic mas- 



ter equation ( 4.10 ) that it takes Gaussian states to Gaus- 
sian states. This will be true for any diffusive stochastic 
master equation that is at most second-order in x or y. 
The sig nifica nce in this case is that we can again use the 
ansatz (3.10), and we need only the equations of motion 
for the five moments. We find the following equations (to 
be interpreted in the Ito sense [p8[) 



dfiiQ/dt 
dfiQi/dt 
dp.2o/dt 



dfio2/dt 



dfiii/dt 



-^iw + Re {Co(0 [M2o - 1 + *Mii] + CUtWl + iy[p2o] + Git) [fm + i]} (4.13) 

-XMio + Re {Coit) [JM02 - i + Mil] + CKtWl + i^im - «] + Q{t) [M02]} (4.14) 

2 - 2^20 - Re [(m20 - 1)^ + /i?i + (1 + t^)/i20 + M?i + 1 

+wSo(M20 - 1 + «Aiii)^ + '"11(1 + l^)^J.lo + ^22(^11 + «)^ 

+ 2uoi\/l + V (m20 - 1 + «/^ii)/^20 + 2^02(^20 - 1 + i/^ii)(/^ii +i) + 2u*2\/l + i^(/iii +^)^J.2o] /2 (4.15) 

-2xMii + 2 + i^ - Re [(/X02 - 1)^ + Mn + (1 + '^)(m?i + 1) + M02 

+Uoo{ipo2 - i + A^ll)^ + Wll(l + '^)(Aill - if + "22/^02 

-I- 2uqiV1 -I- V (i^02 - i + ^ill){^ill - i) + 2uS2(*M02 - i + Mii)/-i02 + 2^*2^1 + v {pn - i)/io2] /2 (4.16) 

-Mil - XM20 - Re {(m20 - 1 + ipii){-ipo2 + i + Mn) + (1 + ^)(Mii - «)M20 + Mo2(mii - i) 

+Uoo(M20 - 1 + iMii)(iM02 - i + Mil) + Mn(l + i^)M2o(mii - «) + '"22Mo2(mii + 

-I- UqiVI + v [(m20 - 1 + iMii)(Mii - «) + M2o(«Mo2 - « + Mn)] + "12VI + t^ [M20M02 + (mh + «)(mii - «)] 

+ w52[(«M02 - i + Mii)(Mii + «) + (M20 - 1 + «Mii)M02]} /2 (4.17) 



C. The Stationary Solutions 

From these equations we see that the evolution of the 
second order moments /i2o, M02, Mn is deterministic. This 
means that for a given unraveling hi the stationary en- 
semble will consist of Gaussian pure states all having 
the same second order moments. They are distinguished 
only by their first order moments x — /iio, y — Moi; which 
therefore take the role of the index n in Eq. (2.11). The 
different ensembles themselves are indexed by another 
pair of n umbers, Mii:M20j which play the role oi lA in 
Eq. ( 2.11 ) . We do not need /io2 because the purity of the 
unraveled states implies that 



However, it should be noted that the mapping from hi to 
Mil: M20 is in general many-to-one as discussed below. 

We now introduce a new notation for the second order 
moments, 



a = M02 ; /3 = Mil ; 7 = M20, 



(4.19) 



The different ensembles are now indexed by the pair /?, 7. 
Of course not all pairs (3, 7 correspond to physically real- 
izable ensembles. Since the ensemble we are considering 
has evolved to a steady state ati = 0, the only valid 



2 1 

M20M02 — Mil — ^■ 



pairs must satisfy Eqs. (4.15)^(4.17) with the left-hand 
sides set to zero. This gives three simultaneous equa- 
tions that, on splitting Uij into real r^j and imaginary 



(4.18) hij components, can be written as 



1 _ ^ _ (1 + j,/2)-y2 _ ^2 ^ ^^^[(^ _ ^Y - l3^]/2 + rn(l + ly)-^^ /2 + 7^22 (/J^ - l)/2 

+ /Joo/3(7 - 1) + h22p 



-2xP+{l + v/2){l-P^) 



+ 7'oi\/l + 1^7(7- 1) +^'02(7- 2)/3 + ri2\/l + 1/7/? 

+ /ioi\/TT^7/3 + ^02(^^ + 7- l) + ^i2\/l + J^7 
^2 + a = roo[/32 - (a - l)2]/2 + rn(l + v){P^ - l)/2 + r22a' /2 

+ /ioo/3(a-l) + /jii(l + J^)(-/3) 

+ roi%/TT^(/32 + a - 1) + ro2/3a + rl2^/l + v|3a 



(4.20) 



(4.21) 



+ /ioi\/rT^(a-2)/3 + /io2(a- l)a + hi2\/l + v {-a) 
-XI - a/3 - (1 + J^/2)7/3 = roo/3(7 - a)/2 + rii(l + i/)7/3/2 + r22/3a/2 

+ /ioo[/3' + (a - 1)(7 - l)]/2 + /iii(l + i^)(-7)/2 + /i22a/2 

+ roW\ + v-il3 + ro2 [/3^ + 1 + (7 - 2)a]/2 + +ri2^Jl + v (a-f + /J^ + l)/2 

+ /loi VTT^ [/32 + 1 + (a - 2)7]/2 + /io2/3a, 



(4.22) 



where a is to be read as (1 + 0^)1^- 

These three equations are nonUnear in /3, 7 but hnear 
in the 12 real variables rij,hij. This means that if the 
values of 7 and (3 are known then the three equations can 
be solved for r-y , hij . Since there are only three equations 
for the 12 unknown variables, the resulting linear system 
is non-singular and an (uncountably) infinite number of 
solutions are possible. We denote the family of such so- 
lutions F^ = {r'lf,h'^^ ■■ ij = 0,1,2}, indexed by C- 
Physically this arises because many different unravelings 
U may lead to the same steady state ensemble /3, 7. The 
question of whether a given pair of values of 7 and /3 
represents a physically realizable state then becomes the 
problem of determining whether any of the solutions F^ 

for the correlation coefhcients uj; = r-J -|- ih-J satisfy 

''J ''J ''J 

the condition ||u(«)|| < 1 in Eq. ( |2.18D . This problem 
can be solved by finding the solution Fs that gives the 
smallest value for ||u||, and checking if this is less than 
one. 

The above method determines the boundary between 
those ensembles that are physically realizable and those 
that are not by finding, explicitly, the parameters of the 
unravelings that satisfy ||u*-^'|| = 1. There is an alter- 
nate, but equivalent, approach Q based on the central 
idea of Ref. |Q. This allows one to take an arbitrary 
ensemble and check whether it is possible, by monitoring 
the environment, for the state of the system to be re- 
stricted to members of the ensemble over arbitrary time 
intervals. The ensemble is physically realizable if, and 
only if, this can be done without changing the ensemble 
average dynamics. The advantage of this alternate ap- 
proach is that the parameters of the unraveling need not 
be calculated explicitly and so the computational task 
can be greatly reduced. Moreover it is possible to find 
the boundary between physically realizable and non re- 
alizable ensembles in a closed analytic form. The details 
are tangential to the scope of the present paper and are 
explored elsewhere [Q. We note here however that the 
PR region is given by /3 and 7 values satisfying 7 > 
and 



We have verified this analytic result with numerical so- 
lutions obtained using the former approach, for all cases 
presented below. 



D. The Stationary Ensemble 

Th e st ationary solution of the linearized master equa- 
tion ( |3.8| ) has a Wigner function which is independent 
of phase (y) and has the following amplitude (x) depen- 
dence: 



Wssix) ex (27r)-i/2gxp(_a;2/2). 



(4.24) 



This is as expected from the stationary solution of the 
full master equation, Eq. ( |3.6[ ). A flat phase distribution 
linearizes into a flat y-distribution. 

A s shown above, the long-time solution of the SME 
(4.9) is an ensemble of Gaussian pure states in which the 
second order moments /i20, A*ii, M02 are identical in all 
members of the ensemble, but x = /iio and y — /ioi are 
allowed to vary. The ensemble is thus represented as 



^'' = {(pL,^l'5):S,ye5R}, 



(4.25) 



where the second order moments of the pure state P!^y 
are determined by the unraveling U. 

The weighting function p^ ^ for the members of the 
ensem ble i s Gauss ian. This follows from the fact that 
Eqs. ( 4.13 ), ( 4.14 ) for x and y describe in steady state 
(where the second-order moments are constant) a two- 
dimensional Ornstcin-Uhlenbeck process |g^. Such a 
process has a stationary probability distribution that is 
Gaussian. 

Rather than deriving this stationary Gaussian distri- 
bution pj y from the Ornstein-Uhlenbeck process we can 
derive it more simply by noting that it must satisfy 



Psi 



= / dxdyp^yP^y. 



(4.26) 



(2x/3 - 2 - y){2 - 27) - (/3 + xif > 0. 



(4.23) 



This is guaranteed by the fact that the SME is equivalent 
to the master equation on average. Evidently y should 
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always have a flat weighting distribution, and x should 
have the weighting distribution 

p"(S) = [2^(1 - Ai2o)]"'/' cxp [-xV2(l - ^l2o)] ■ (4.27) 
This ensures that 

Wssix) ^ j dy j dxp'^{x)Wly{x,y), (4.28) 

where W^y{x,y) is the Wigner function of Pj^y. 

V. PR ENSEMBLES FOR THE (ATOM) LASER 

In this section we present our results for the physically 
realizable ensembles for the (atom) laser. 

A. Realizing the Number State Ensemble 

Before turning to the effect of varying the dynamical 
parameters x ^-nd i> we briefly return to the physical re- 
aliza bility of the number state ensemble. We showed in 
Sec. IV A above that this ensemble can be realized by 
a discontinuous unraveling. The analog of the number 
states in the linearized regime we have been consider- 
ing are the inflnitely squeezed states with 7 = /3 = 0, 
a = 00. We expect that these states should be PR using 
a CM unraveling. This expectation is m et, in t hat t hese 
state parameters are a solution of Eqs. ( 4.2C| )-( 4.22 ) for 
Moo = Ij Wii = 1, U22 = —1, and all other u^j, — 0. 



B. Varying x ^vith v — 

First we present the results showing the effect of vary- 
ing X for fixed v = 0. As we have established above, a 
PR ensemble from a CMU can be represented by the pair 
of numbers 7, (3. Thus the set of all PR ensembles can be 
represented by a region in 7 — /? space [0,1] x (—00 , cxd). 
The boundaries of this region, given by Eq. (4.23), are 
shown in Fig. m for various values oi x- A number of 
features of this plot are evident. First, for any non-zero 
value of Xi the coherent state ensemble is not PR. Sec- 
ond, as X increases the PR ensembles become increas- 
ingly removed from the coherent state ensemble. Third, 
the boundary of the PR ensembles is asymmetric in (3 for 
X > 0, with a larger negative f3 region. 

The first point can easily be proven analytically. Co- 
here nt sta tes are given by a = 7 = 1 and /3 = for which 
Eq. ( 4.23 ) gives — x^ > 0. That is, coherent states are 
physically realizable only for x = 0. 

We quantify the second point by defining the closest- 
to-coherent (CC) ensemble as that for which the states 
have maximum overlap with a coherent state. The over- 
lap of two Gaussian states with the same mean ampli- 
tudes and covariance parameters a,/3,j and ai,/3i,7i is 



2/v/(ai+«)(7i+7)-(/3i+/3)2 



(5.1) 



If one of these is a coherent state, with ai =71 = 1, 
Pi = 0, this reduces to 



/2-|-a + 7. 



(5.2) 



Thus, to find the closest-to-coherent ensemble we simply 
find the minimum a +7 = 7+(H-/3^)/7 in the PR region 
of 7 — /? space. 

The closest-to-coherent ensemble for each value of x is 
represented in Fig. H^ as a filled circle on the boundary of 
the respective PR region. The states in these ensembles 
become more squeezed (7 —^ 0) and have a greater x — y 
covariance as x increases. This trend is shown in more 
detail in Fig. where we plot the parameters a, (3 and 
7 for the closest-to-coherent PRE as a function of x- By 
finding the min imum of 7 -|- (1 -|- /3^)/7 subject to the 
constraint Eq. (4.23) and expanding about 7 = and 
1/x = we find the parameters of the CC ensemble for 
large x scale as 



CC _ ^^1/2 



CC 



33/4- 

2 



X 



~ 31/4 
iCC ^ 1 

- 31/2 



-1/2 



(5.3) 
(5.4) 

(5.5) 



Also plotted in the figure are two lines representing x^^^ 
and x^^^^ foi' comparison. One can clearly see the 1/2 
power law scaling for a and 7. 

The third point, i.e. the increasing asymmetry of the 
PR regions in Fig. nl is due to the self-energy of the 
condensate embodied by the term containing (a^a)'^ in 



Eq. (3.4). In the Wigner phase-space representation, this 



term by itself produces a 'phase shearing'; that is, the an- 
gular velocity of the point {x,y) depends on the distance 
~ y/JI{l-\-x) from the origin [Bq|. In our linearized model 
of the atom laser, the effect of this term is to shear t he cir - 
cular contours of a coherent state into ellipses. Eq. ( 3.19| ) 
indicates that these ellipses have a negative covariance. 
When monitoring the reservoirs it will therefore be eas- 
ier to realize states with a negative covariance. Hence, 
the PR regions become more asymmetric allowing more 
negative-/? regions as the nonlinearity parameter x in- 
creases. 



C. The effect of nonzero v 

Non-zero values of v, as defined in Eq. ( ^.9[) , corre- 
spond to the presence of excess phase diffusion, which 
will tend to overcome the phase-shearing effect. This 
makes it easier to physically realize states that are closer 
to coherent states. In Fig. |^ we plot the boundaries of 
the PR ensembles for 1/ = 10 for the same set of values 
of X as in Fig. |l|. The CC ensembles are also shown as 
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filled circles. The PR regions are generally broader as 
expected, and this allows the CC ensemble to be closer 
to a coherent state than for corresponding x values in 
Fig. I 

The parameters for the CC ensembles for v = 100 as 
a function of x are plotted in Fig. |j. Comparing with 
Fig. H we note that the presence of the excess phase dif- 
fusion in Fig. y allows ensembles very close to coherent 
states (i.e. with a«7~l,/3~0) for x up to of order 
z/^/^. This can be verified analytically from Eq. ( 4.23 ). 
However, as the value of x increases beyond this to of 
order i/, the effect of the non-zero i^ value becomes less 
significant and the curves approach the same asymptotes 
as in Fig. ||. 

The physically-realizable region for x — ^ includes the 
point /? = 0, 7 = 1 for all values of i/. Hence the closest- 
to-coherent PR ensemble is trivially an ensemble of co- 
herent states in this case. The situation is different for 
nonzero x- Fig- B shows the parameters for the closest- 
to-coherent PR ensemble as a function of u for x = 100. 
For v Ki the values of a, (3 and 7 are approximately the 
same as the corresponding values at x = 100 in Fig. 0. 
However, as v increases much larger than %, the effect 
of the self-energy term become less significant and the 
phase diffusion begins to dominate. Then, for 1^ ^ x^ 
the closest-to-coherent ensemble approaches a set of co- 
herent states as a, 7 ^ 1. 



D. Comparison with quantum state diffusion 

The unraveling given by quantum state diffusion 
(QSD) is more restrictive than that of the general contin- 
uous Markovian unraveling treated here. Specifica lly, fo r 
the QSD unraveHng, a, /3 and 7 must satisfy Eqs. (4.20)- 



(4.22) for 



^ij 



'ij 1 tfiij 



instead of any Uij fulfilling 
||u|| < 1. We find this yields the analytic solutions for 
the QSD ensemble 



(5.6) 
(5.7) 
(5.8) 



QSD 


1 + Vl-8x/3 + 4M(l-/32) 




2 


QSD 


-l + v/l + 4M(l-/32) 




2M 


QSD _ 


{-l + UI-F)x + VG-E 



4(x2 + M) 



where 



M = 1 + iy/2 
E = (24Af - 2)x^ 



32M-' + 8M^ 



i^ = 4V(Af + 1/4)2 + x^ 
G= 2{4:M^ + x^)F . 



(5.9) 
(5.10) 
(5.11) 
(5.12) 



The crosses in Figs. W and represent the QSD ensembles 
for the same set of x and ly values as the CC PR ensem- 
bles. The corresponding value of x for the crosses reduces 
from left to right. One immediately notices that the QSD 



ensembles lie well inside of the PR boundary indicating 
that, for moderate x and v values, the QSD unraveling is 
significantly more restrictive than the general continuous 
Markovian unraveling explored here. Moreover, the QSD 
ensembles are more squeezed (smaller 7 values) than the 
corresponding CC ensembles. 

We note that the QSD ensemble is significantly 
squeezed even for the ideal photon laser limit of x = 
1/ = for which the QSD ensemble is given by a = 
{V5 + l)/2 « 1.62, /3 = and 7 = (v^ - l)/2 « 0.62. 
We can trace the origin of this squeez ing as follows. The 
second term on the right side of Eq. ( |3.4| ) represents the 
output coupling of the laser. As mentioned above, QSD 
corresponds to equal-efficiency homodyne detection of a 
pair of orthogonal quadratures. Thus, in QSD the moni- 
toring of the output will tend to localize the state of the 
laser onto a coherent state. No squeezing can therefore 
originate from this term. The squeezing must therefore 
originate from the nonlinear amplification process repre- 
sented by the first term on the right side of Eq. (pA). 
Indeed, the nonlinear amplification restricts the ampli- 
tude noise through depletion of the source. In our lin- 
earized model, this corresponds to restricted noise in x. 
Evidently, the monitoring of the reservoir modes associ- 
ated with the amplification is a partial measurement of 
X and this leads to the squeezing of x. 

It is interesting to compare this with the general con- 
tinuous Markovian unraveling (CMU) treated in the pre- 
vious subsection. This is less restrictive than QSD since, 
for example, it allows the unbalanced monitoring of two 
quadratures of the output field. In particular, a correla- 
tion value of Woo = ~1 corresponds to the monitoring of 
just the y quadrature. This would tend to localize the 
state of the laser mode onto a state with reduced y fluc- 
tuations and thus counteract the x-quadrature squeezing 
effect from the nonlinear amplification. Similar remarks 
apply to unraveling the gain process itself. The net effect 
is that the general continuous Markovian unravelings can 
physically realize coherent states for x = ^^ = whereas 
QSD does not. 

Despite these differences, the a and 7 scaling laws for 
the QSD ensemble follow the same X^^''^ power laws as 
the closest-to-coherent ensemble although with a differ- 
ent prefactor. In Fig. O we plot the parameters for the 
QSD ensemble for 1^ = as a function of x- Compar- 
ing with Fig. 3 we note that the QSD ensemble begins 
more squeezed for small Xj but for large x the two en- 
sembles ap pro ac h si milar degrees of squeezing. In fact, 
from Eqs. (5.(:)-(5.S) we find the scaling laws 



^QSD ^ _^ 



(5.13) 
(5.14) 
(5.15) 



which should be compared with Eqs. (5.3)-(p.5|). 

In Fig. we plot the parameters for QSD ensemble for 
X = as a function of v. The QSD ensembles are highly 
squeezed for increasing v and, indeed, we find 
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aQS^c 


\/2 


^QSD^ 


. V2^-i/2 


^QSD ^ 


= . 



(5.16) 

(5.17) 
(5.18) 



This is perhaps surprising given that one does not usuaUy 
associate enhanced squeezing with large phase diffusion. 
However, the monitoring of the reservoir corresponding 
to the phase diffusion is effectively an incomplete mea- 
surement of the variable a^a, which, in our linearized 



model Eq. (3.8), is represented by the term vV\xl'2\. 
The monitoring therefore tends to localize the state of 
the laser onto an eigenstate of x. The strength or rate 
of these measurements increases with v. In QSD there is 
no mechanism to counteract the associated squeezing of 
the a;-quadrature, and so the squeezing increases with v. 
In contrast, the general continuous Markovian unraveling 
allows unbalanced monitoring of all baths. In particular, 
with u\\ = —1, the phase diffusion is unraveled as a pure 
noise process (stochastically changing the phase of the 
state, but yielding no information about it). This al- 
lows the closest-to-coherent CMU ensemble to comprise 
of coherent states for the same parameters as for Fig. 0. 



VI. DISCUSSION 
A. Summary 

The atom laser, even under with the simplifying ap- 
proximations we have made, is an open quantum system 
with rich dynamics. Some aspects of the dynamics, such 
as excess phase diffusion (parametrized by v) and phase 
dispersion caused by atomic interactions (parametrized 
by x)j do not affect the stationary state. That is because 
the stationary state is a Poissonian mixture of number 
states. In this paper we have investigated the represen- 
tations of this mixed state as ensembles of pure states. 
The diagonal representation (number states) is one such 
ensemble, and the random-phase coherent state ensem- 
ble is another. Although mathematically equivalent we 
have found that such representations are not physically 
equivalent, as only some of them can be physically re- 
alized through monitoring the system. Moreover, the 
dynamical parameter x, which does not affect the sta- 
tionary state at all, radically affects which pure state en- 
sembles are physically realizable (PR). In particular, for 
any x t^ Oj the ensemble of coherent states with unknown 
phase is not PR. 

As the nonlinearity x is increased, the PR ensembles 
become increasingly removed from the coherent state 
ensemble. To be specific, the ensemble of states that 
are closest to coherent states consists of states that are 
amplitude squeezed (but slightly rotated), with a phase 
quadrature variance increasing as 



As X increases the closest-to-coherent (CC) ensemble be- 
comes more squeezed until eventually the linearization 
leading to the above result breaks down. This indicates 
that is not possible to physically realize an ensemble with 
a well-defined coherent amplitude for a x this large. This 
occurs when cP'^ ~ /i, in other words x ^ M^- Note that 
this is larger than the critical value x ~ \^^'^ at which 
the laser becomes incoherent, according to the analysis 
of Sec. pT^ . 

The situation is quite different in terms of the excess 
phase diffusion parameter v. As v increases (with x = 0) 
the coherent state ensemble remains PR. This is true even 
when V > [1? ^ the value at which the laser becomes inco- 



herent, as shown in Sec. IIIC. Moreover, phase diffusion 



tends to undo the nonlinear effects of the self energy. In 
the limit z. — > cx), the coherent state ensemble is PR for 
any finite value of x- 



B. Interpretation 

In Ref. p|, the coherence condition for a laser, that 
the output flux be much greater than the linewidth, was 
motivated by the requirement that the laser have a well- 
defined phase. This follows from the following argument. 
The laser phase remains fairly constant over the coher- 
ence time (the reciprocal of the linewidth). However this 
phase only has meaning if it can be measured, and this 
requires a macroscopic field (i.e. many bosons) to be 
produced in the ou tput over one coherence time. As de- 
rived in Sec. IIIC, this condition requires x ^ [^ ^^nd 

From the results of this paper there seems to be a prob- 
lem with this motivation for this definition of coherence. 
There are values of x between /z^" and /i^, and v be- 
tween ix^ and oo, for which the atom laser is not coherent 
and yet for which it is possible to physically realize laser 
states with well-defined coherent amplitudes. 

The resolution of this problem is straight-forward for 
the case of large v. The motivation in Ref. [B| relied upon 
a measurement of the phase jrora the laser output. By 
contrast, the ensembles we have considered in this paper 
are physically realized by monitoring all of the reservoirs 
of the laser. In particular, that means monitoring the 
reservoirs that produce the excess phase diffusion v. If 
we only allow for monitoring of the output of the laser, 
the stochastic master equation will not preserve purity. 
After linearization, the following equation results 



dp^ {l/'i)dt{V[x + iy]p- 
+ H[i{xy + yx)/2]^ 
+ {l/2)n[dW*{t){x 



a 



u)V[x\p 

n[~ixx^]]p 

Viy)]p. 



T^[y]p 



(6.2) 



„cc 



X 



1/2 



(6.1) 



Here there is only one stochastic term, from monitoring 
the laser output. The best strategy for trying to realize 
states with well-defined coherent amplitudes is clearly to 



13 



measure the phase quadrature of the output. This corre- 
sponds to dWdW = —dt. 

Under these conditions, the differential equations for 
the second-order moments of the conditioned state are 



fJ'20 

All 



2 - 2/120 
-Mil - XM20 



M?i, 



(M02 - 1)M11, 

M02 = -2xMii + 2 + z/ - {fiQ2 - 1)^- 
If we set X — 0, the steady-state sohitions are 



^J■20 = 1, 
Mil = 0, 
/i02 = 1 + V2" 



(6.3) 
(6.4) 
(6.5) 



(6.6) 
(6.7) 
(6.8) 



In the Hmit of large v (which is the potential problem 
area), the phase quadrature variance scales as v^^^. The 
states lose their coherent amplitude as the linearization 
breaks down at a = /io2 ^ /i- That is to say, at ly ^ fi^. 
This is precisely the regime identified in Sec. Ill C as that 



for which the laser output loses its coherence. 

Unfortunately (or perhaps fortunately from the point 
of view of provoking new concepts), a similar analysis for 
large x does not hold. Instead, with i^ — 0(1) and x ^ 1 
the solutions of Eqs. (6^)-(^^) are 



1/2 



/i20 ^ 2'^\ 

fill ~ -V2, 

M02 ^ 2'/^x'/\ 



(6.9) 
(6.10) 
(6.11) 



This is an extremely sheared state, with phase quadra- 
ture variance scaling as x^^^- K loses its well defined 
phase only for x ~ M^ j which is the same scaling as found 
above when all the reservoirs were unraveled. In partic- 
ular, for /i"^'^ < X < M'^j measuring the output has deter- 
mined the phase of the laser even though this should not 
be possible by the argument in Ref. because the flux 
is less than the linewidth. 

The difference between large v and large x can be un- 
derstood as follows. There are three Lindblad terms in 
the linearized master equation ( |3.8| ). When v — they 
are all of roughly the same size. Thus restricting the mon- 
itoring to just one of the three reservoirs (the first one, 
the output) has relatively little effect on the conditioned 
states. It is much like monitoring all reservoirs, but with 
a reduced efficiency. Indeed, the conditioned state in this 
case is not far from a pure state, with /i2oMo2 — Mii = 2 
(compared to 1 for a pure state). By contrast, with i^ 
large the phase diffusion Lindblad term is much larger 
than the other two. Then if one is only able to moni- 
tor the output one is necessarily losing most of the in- 
formation about the system. This leads to qualitatively 
different conditioned states, with much reduced purity 
(M20M02 - Mil - v^ > !)• 

The existence of the regime n^'^ < x < M^ where the 
laser output is incoherent, but where the phase can in 



fact be determined suggests that the concept of coher- 
ence time is more subtle than the standard definition in 
terms of the first order coherence function used in Ref. ||] 
and in Sec. [II C above. The coherence time is also used 



to define whether or not the laser beam is Bose degen- 
erate, and, as discussed in Ref. [g|, the criterion is the 
same. That is, the output is Bose degenerate if and only 
if many bosons come out "with the same phase" (that is, 
within one coherence time). Thus the present paradox 
has implications that go beyond the present discussion, 
and impact on concepts like Bose degeneracy as well, as 
will be discussed below. 



C. Conditional Coherence and Conditional 
Degeneracy 

One way to understand the above results is that the 
atom laser for fr^^^ < x < M^ is "conditionally coher- 
ent". The standard coherence condition x < M^^'^ can 
be derived from the requirement that ((J0)^(i)) < 1 at 
t — l//i, the time between atoms in the output. Here 
((5(/i)^(i) ~ 2/^(i)/4/i is the phase variance of the state at 
time t, which was a coherent state at t = 0. That this im- 
plies the condition x < M"^^^ can be seen simply as follows. 
For X large and for a time as short as l//x, the irreversible 
evolution can be ignored and the phase uncertainty is due 
to the Ca^'a^aa Hamiltonian. For the linearized theory, 
this turns into the Hamiltonian x(2^/2)^, where x is the 
amplitude quadrature. This causes the phase quadrature 
to change as 



2/(t)=y(0)-xte(0), 



(6.12) 



where the mean frequency shift has been removed as has 
been consistently done before. For a coherent state of 
zero mean phase we have y{0) = 0, {y{0)^} = 1 and 
x{0) = 0, {x{Of) = 1. Thus for t = l//i we get 



{[yit)-mf)^i + x't'^i + xVti 



(6.13) 



This is of order 4/i (indicating the loss of coherence) for 

X-M^/^. 

The coherent state is the most convenient state to use 
for this calculation, as explained in Sec. HI C . But of 



course it is also possible to represent the atom laser as 
a mixture of states with smaller amplitude uncertainty 
than a coherent state, and, as we have seen, to physi- 
cally realize such ensembles. The average result must be 
the same, but the details are different. Consider a mini- 
mum uncertainty pure state with V = (^[x{0) — x{0)]'^) = 
1/ (j/(0)'^), where the initial mean phase has again been 
taken to be zero. The mean phase evolves as 

m = -xtxiO), (6.14) 

and the phase quadrature variance as 

{[yit)-m?) = yV + x't'V (6.15) 
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To reproduce the stationary state which has a unit vari- 
ance, we must consider an ensemble of different vahies 
for 5(0), with mean zero and variance 1 — V. Thus the 
total phase variance over the ensemble, 

{[y{t) - m?) + E [f] - i/v + x^ey + x'^'e[s(o)2] 

= l/F + x'i' (6.16) 

cannot be less than that from a coherent state (with 

In this picture, the increase in the phase uncertainty 
is the sum of an intrinsic phase uncertainty increase and 
that due to an uncertainty in the frequency of the field. 
The former is due to an initial quantum uncertainty V 
in the amplitude quadrature, and the latter to a classical 
uncertainty 1 — V^ in the initial mean amplitude quadra- 
ture. The loss of coherence is thus partly due to the addi- 
tion of different interference terms oscillating at different 
frequencies. For example, interfering parts of the output 
field separated in time by t would give a different inter- 
ference pattern depending on the frequency. Over a time 
of order unity (the bare decay time) , the mean amplitude 
will sample all possible values so the frequency will also 
vary. The average interference pattern measured over a 
time long compared to this will thus be washed out due 
to the different frequencies, and the experimenter would 
conclude that the output was incoherent if -y^t^ ~ [i for 

t - 1/m- 

If, however, one knows (as the experimenter) the initial 
mean amplitude a;(0), then one knows what frequency to 
expect in one's interference pattern. Then rather than 
simply averaging the interference patterns over some long 
time, one could correct for the mean frequency shift be- 
fore doing the average. Then the only contribution to 
the visibility of the interference patter will be the intrin- 
sic phase quadrature variance 



{[yit)-m?)^yv + x't'V. 



(6.17) 



From this conditional point of view, the laser output will 
cease to be coherent only when 



Solving for x gives 

x-^iVv-H^^i-v-^). 



(6.18) 



(6.19) 



To maintain coherence for the largest possible x, we min- 
imize this with respect to V to get 



X^2/i2 



(6.20) 



at F '^ 1/2/i. This is the upper limit of the region 
^3/2 ^ ^ ^ ^2 .^^]^g]-g a^ well-defined coherent amplitude is 
physically realizable but the output is not coherent in the 
usual sense. Now we can see that a physical realization 
giving the well-defined coherent amplitude in this regime 
(such as that giving the closest-to-coherent ensemble) is 



precisely what is required to recover coherence, in a con- 
ditional sense. 

The concept of conditionally coherent goes hand-in- 
hand with that of conditionally Bose degenerate. Un- 
der the standard definition, the atom laser output in the 
regime //■^/^ < x < M^ is not Bose degenerate. Specifi- 
cally, there is no mode that can be identified a priori in 
the output and that has a large mean occupation number. 
But under an unraveling of the atom laser dynamics, such 
a mode can be identified in this regime: it is a mode cor- 
responding to the frequency which can be inferred from 
the knowledge of the amplitude of the condensate. As 
with the case of conditional coherence, a new mode will 
have to be chosen after a short time, since the frequency 
explores the full range on a time scale of order unity. 
But at a particular instant of time, the knowledge ob- 
tained from monitoring the reservoirs of the system (or 
even just the output, as seen above) is sufficient to allow 
a highly-occupied mode to be identified. 

We can perhaps clarify the concept of conditional Bose 
degeneracy as follows. Consider a system with TV modes, 
and N particles. The multiparticle state 



p- lA^l, 02,03,..., Oat) (iVi, 02,03 



,0 



ATI 



(6.21) 



is clearly Bose-degenerate, just as the state 

p=|li,l2,l3,...,lAr)(li,l2,l3,...,lw| (6.22) 



is not. But what about the state 

N 



P 



— N 2_^ I- ■ • Om-2, 0„i_i, Nm, Om+1, 0„i+2 ■ 



m—1 



X (. . . 0™_2, 0™_i, TV™, 0™+i, 0™+2 . . .| ? (6.23) 



The mean occupation number of any mode is clearly one, 
so it is not Bose degenerate in the usual sense. But also 
clearly if one had access to this state then after finding 
a single particle, one would know in what state the re- 
maining A^ — 1 particles would lie. Thus the state would 
have become conditionally Bose degenerate. We believe 
that the above state is a good toy description of a short 
section of the output of an atom laser in the interest- 
ing regime of /i'^" < x < M^j where the different modes 
represent different frequencies. 

Finally, it is interesting to note that by employing feed- 
back based on QND atom number measurements, it is 
possible (within the current atom laser model) greatly to 
reduce the linewidth [Q . Specifically, the linewidth may 
be reduced by a factor of order /i^/^, and the coherence 
(in the conventional sense) of the laser extended from 
X ^ M^^^ to X ^ A^^. This is not quite an exact paral- 
lel with the above results, because the feedback is based 
on a measurement that adds extra phase diffusion (v) 
term, that is not required in the above analysis. (This 
QND measurement is introduced because it is a num- 
ber measurement, and so is more easily realized than the 
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phase-sensitive measurement necessary in the above anal- 
ysis.) Nevertheless, it still illustrates the general princi- 
ple stated in Ref. |Q, that "the practical significance of 
[conditional analyses] is that conditioning is realized by 
feedback." 



D. Experimental Implications 

It is clear that many interesting questions relating to 
the coherence of an atom laser, the physical realizability 
of a coherent state ensemble, the coherence of the out- 
put, and the conditional coherence of the output, depend 
upon the value of %. This prompts the question: what 
value has this parameter in experimental atom lasers? As 
discussed in the introduction, a number of experimental 
groups have realized Bose- Einstein condensates with out- 
put coupling jT^-pQl- A CW atom laser would have to 
incorporate a mechanism for replenishing the condensate 
so that the output coupling could continue indefinitely. 
Nevertheless we can take these experiments as a possible 
indication for the parameter regime in which an atom 
laser may work. The figures below are derived by setting 
the bare linewidth k of the laser equal to the reciprocal 
of the lifetime of the condensates in the experiment, and 
the mean atom number fi equal to the initial occupation 
number of the condensate. The excess phase diffusi on t / 
we have i gnor ed, and we have calculated x using Eq. (3.3) 
and Eq. (|3j|). 

Most current experiments work in the regime where 
the ratio of the kinetic energy to the interaction energy 
is very small W8[: 



QAn'^muj^'^al 



2/5 



<1. 



(6.24) 



Here m is the atomic mass, to is the mean trap frequency, 
and Os is the scattering length as in Eq. (3.3). In this 
regime the Thomas-Fermi approximation can be made, 
allowing us to evaluate x analytically as 



4 f225fi'^muj^al 



7k \ 



1/5 



(6.25) 



The values of x using the parameters of three recent ex- 
periments are compared in Table I. The MIT experiment 
||17| represents the first "pulsed atom laser" , a quasicon- 
tinuous output coupling [ 19| was demonstrated at NIST, 
and the MPQ experiment |[20[| demonstrated a continuous 
output coupling. 





MIT 


MPQ 


NIST 


Proposed 


X 


910 


1800 


50 


990 


IT 


4.1 X 10^ 


2.1 X 10^ 


5.7 X 10^ 


2.9 X 10^ 


i/e 


6.0 X 10^ 


4.0 X 10^ 


8.0 X 10^ 


2.0 X 10^ 


K/e 


12 


0.57 


810 


2.0 


OJmin/K 


1.1 


4.8 


0.8 


22 


UJmin/i 


14 


2.7 


640 


44 



Table I. Parameters for recent (and proposed) atom laser 
experiments at various institutions 

All X values are in the x ^ 1 regime on which we have 
concentrated in this paper. Thus if these experiments 
could be run with the same output coupling but with 
continuous replenishment of the condensate, the closest- 
to-coherent ensemble that could be physically realized 
would be highly amplitude squeezed. From Eq. (p]j), 
with X = 1000 the standard deviation of the amplitude- 
quadrature of these states would be about 0.2, compared 
to 1 for coherent states. Thus it seems that it is wrong 
to think of an atom laser as being in a coherent state. 

Despite the banishing of the coherent state descrip- 
tion, truly continuous versions of the experiments ana- 
lyzed above would produce an unconditionally coherent 
(Bose degenerate) output. That is because the calcu- 
lated values of X are always much less than /i'^/^, so that 
Eq. ( |3.30|) above is satisfied. Interestingly, we can recast 
this condition in terms of the output flux I ^ Kfi (atoms 
per unit time) as 



/> 1.61w 



1/11 



(6.26) 



This inequality depends very weakly on the dimension- 
less quantity in brackets because of the 11th root. For 
the above three experiments this 11th root averages to 
0.16, and ranges only from 0.13 to 0.21. Hence we can 
state the coherence condition for an atom laser in terms 
essentially independent of the species and decay time as 
/ > 0.26a; or 



/>T- 



(6.27) 



That is, there should be many atoms emitted into the 
laser beam per oscillation period T = 2tt/u; of the trap. 
This is such a simple rule of thumb that it should be 
useful, but it must be remembered that there is no di- 
rect physical connection between the flux and the trap 
frequency. This result is simply a numerical coincidence 
arising from the various physical parameters for atomic 
Bose-Einstein condensation in typical traps. The second 
row of the table shows that this condition is clearly satis- 
fied for the parameters of the three experiments and this 
suggest that the output field of our model atom laser 
would be degenerate. 

The actual degree of degeneracy D of the output field, 
that is the number of atoms per output frequency mode, 
is given by the quotient I /£ of the flux / and the linewidth 
£. The linewidth for the atom laser model we are consid- 
ering is given in Ref. B^ as 



2kx/\/2¥JI 



for X < V^fi/ir 
for X > y/siJJiT 



(6.28) 



The third row of the table shows that, for the same pa- 
rameters as the experiments, the output field of the atom 
laser model is highly degenerate. 
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It is interesting to compare the linewidth of the output 
field i with the bare cavity linewidth n. The action of 
the pump tends to reduce the linewidth far below k in 
the s ame manner of an optical laser [the x ~^ limit of 
Eq. (6.28)]. In an atom laser, however, the nonlinear- 
ity converts intensity fluctuations into phase fluctuations 
and thi s tend s to broaden the linewidth [the x ~^ oo limit 
of Eq. ( |6.2§| )]. Table I shows a range of values of the ra- 
tio K,/i from below unity (line broadening) to well above 
unity (line narrowing) for the parameters of the experi- 
ments. We can write k — I / jjl and t = I /D and so the 
ratio k/^ = -D//J. is also the ratio of the number of atoms 
per output frequency mode to the steady state popula- 
tion in the cavity. Significant line narrowing therefore 
leads to Z? 3> /i, that is, many more atoms per output 
mode than in the condensate. 

Our analysis assumes that we can treat the atomic con- 
densate as a single atomic field mode. We now show 
how this assumption can be justified with realistic ex- 
perimental conditions. Only a single mode is needed if 
the condensate is, at most, only weakly coupled to the 
quasiparticle modes. There are two important ways in 
which this coupling can arise. One is due to the fact 
that the spatial form of the quasiparticle modes depends 
on the number of atoms in the condensate and so fluc- 
tuation in the condensate number will cause an overlap 
between condensate and quasiparticle modes. However, 
provided the fluctuations in the condensate atom number 
occur on a time scale much longer that the dynamics of 
the condensate and quasiparticle modes, the system will 
evolve adiabatically and remain in the condensate mode. 
Thus the first requirement for minimal coupling to the 
quasiparticle modes is 



the table shows the values for a Sodium atom laser in a 
symmetric trap with frequency uj = Wniin = 27r x 25IIz, 
output coupling rate k — 7s^^ and mean atom number 



Wmin/K > 1 



(6.29) 



where Lo-min is the lowest of the trap frequencies. The 
other coupling mechanism is due to the linewidth of the 
condensate mode. In order to avoid adiabatic exchange 
of atoms between condensate mode and quasiparticle 
modes, we need the linewidth to be much smaller than 
the spacing between the condensate mode and first ex- 
cited mode. This difference is simply the lowest trap 
frequency Wmin [Q • Hence the second requirement for a 
single mode analysis is 



> 1 . 



(6.30) 



We have tabulated figures for these parameters in Table I 
for the three experiments and included further data for a 
proposed experiment. The three experiments are clearly 
not operating in the single mode regime as lo^i^/ n or 
'^min/^ or are order unity. So besides not being con- 
tinuously pumped, the experiments also do not satisfy 
the single mode criteria of our model and thus require 
a pulsed, multimode analysis such as that of Ref. [ pO[ . 
However it would not be difficult to achieve single mode 
operation by selecting different, but experimentally rea- 
sonable, parameters. For example, the last column in 



of /i = 10^. Both conditions Eq. ( |6.29D and Eq. ( |6.30| ) 
are satisfied and so the coupling would be minimal in this 



E. Closing Remarks 

It is fitting to end by referring to the very beginning, 
that is, the title of our paper. What does the physical 
realizability of ensembles of pure states say about atom 
lasers, coherent states and coherence? 

First, they estabhsh a basis on which it is possible to 
objectively discuss the existence of coherent states as the 
state for an atom laser. 

Second, they show that these coherent states can only 
exist (that is, be physically realized) for x = (that is, 
in the total absence of interactions between the atoms). 

Third, the existence of pure states close to coherent 
states requires x *^ Ij which is a much stronger condi- 
tion than the x ^ /^'^^^ needed for the laser output to be 
coherent (Bose degenerate). 

Fourth, the existence of states with well-defined coher- 
ent amplitude (that is, with phase variance small com- 
pared to unity) requires x *^ M^i a- far weaker condition 
that that needed for realizing coherent states, and also 
weaker than that required for output coherence. 

Fifth, in the regime /z^^^ ^ X "^ M^j ^ n^w concept of 
coherence (and Bose degeneracy) pertains, that of con- 
ditional coherence (or conditional Bose degeneracy). In 
this regime, knowing which member of physically realiz- 
able ensemble one has at a given point in time allows the 
coherence to be demonstrated, where it could not be in 
the absence of that knowledge. 

Sixth, unlike X) excess phase diffusion u does not de- 
stroy the physical realizability of the coherent state en- 
semble (for X = 0)i and in fact makes it easier to ap- 
proach this ensemble for finite x- 

Seventh, the existence of a regime {v > /i^) in which 
the laser output is incoherent but an ensemble of states 
with well defined coherent amplitudes (indeed, coherent 
states) is physical realizable, does not require a new con- 
cept of coherence. Rather, by restricting the measure- 
ment of the atom laser to the monitoring of its output 
beam itself, the physical realizability of such an ensemble 
is restricted to the coherent-output regime v <$^ fi". 
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FIG. 1. Representation of physically-realizable ensem- 
bles, arising from general continuous Markovian unravelings 
(CMU), for u — and various values of X- The shaded re- 
gions represent values of 7 and l3 (and thus a = (1 -I- 0^)M) 
that can be realized by monitoring. The progressively-darker 
shaded regions correspond to values of x of 0, 1, 4, 16 and 
1000 and are bounded by solid, dashed, dash-doted, dotted 
and dash-dot-doted curves, respectively. The 7, 13 value of the 
closest-to-coherent (CC) ensemble in each region is marked as 
a filled circle on the boundary. The crosses mark the 7, j3 val- 
ues of the quantum-state diffusion (QSD) ensembles for the 
same set of x and v values, with the x values reducing from 
left to right. 
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FIG. 2. The parameters for the physically-realizable en- 
semble that is closest to a coherent ensemble (CC) as a 
function of x with v = 0. The ensembles arise from gen- 
eral continuous Markovian unravelings. These parameters are 
the phase quadrature variance a'^'^ (dotted line), the ampli- 
tude-quadrature variance 7 (dashed line) and the covari- 
ance /3 (dash-dot line) for the members of this ensemble. 
Also shown for comparison are thin solid curves representing 
X'/' and x-'^'- 
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FIG. 4. The parameters for the closest-to-coherent physi- 
cally realizable (CC PR) ensemble as a function of x similar 
to Fig. 2 but here with v = 100. The excess phase diffusion 
allows the realization of states very close to coherent states 
until X ~ ^^''^ 
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FIG. 5. The parameters of the closest-to-coherent physi- 
cally-realizable (CC PR) ensemble similar to Fig. 2 but here 
as a function of v and with x = 100. 
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similar to Fig. 1 but for i/ — 10. 
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FIG. 6. The parameters of the ensemble arising from quan- 
tum state diffusion (QSD) as a function of x with v — 0. The 
labeling follows Fig. 2. 
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FIG. 7. Similar to Fig. 6 but as a function of v with x = 0. 
The thin sohd curves represent values of v^''^ and v^^'"^ . 
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